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ABSTRACT 

We present high-resolution VLA observations of the molecular gas in the host galaxy of the 
highest redshift quasar currently known, SDSS J1148+5251 (z=6.42). Our VLA data of the 
CO(3~2) emission have a maximum resolution of 0.17" x 0.13" (<1 kpc), and enable us to resolve 
the molecular gas emission both spatially and in velocity. The molecular gas in J1148+5251 is 
extended to a radius of 2.5kpc, and the central region shows 2 peaks, separated by 0.3" (1.7kpc). 
These peaks account for about half of the total emission, while the remainder is more extended. 
Each of these unresolved peaks contains a molecular gas mass of ~ 5 x 10 9 M Q (similar to the total 
mass found in nearby ULIRGS) and has an intrinsic brightness temperature of ~ 35 K (averaged 
over the lkpc-sized beam), comparable to what is found in nearby starburst centers. Assuming 
that the molecular gas is gravitationally bound, we estimate a dynamical mass of ~ 4.5 x 10 10 M Q 
within a radius of 2.5 kpc (~ 5.5 x 10 10 M© if corrected for a derived inclination of i ~ 65°). This 
dynamical mass estimate leaves little room for matter other than the detected molecular gas, and 
in particular the data are inconsistent with a ~ 10 12 M Q stellar bulge which would be predicted 
based on the Mbh — ftmigc relation. This finding may indicate that black holes form prior to 
the assembly of the stellar bulges and that the dark matter halos are less massive than predicted 
based on the black hole/bulge mass relationship. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: ISM — quasars: general 
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1. Introduction 

More and more objects have been found at the 
highest rcdshifts in recent years, back to the 'Dark 
Ages' (the epoch of formation of the first lumi- 
nous structures), at z>6 when the universe was 
less than a gigayear old (e.g., Fan et al. 2002, 2003, 
2004, Hu et al. 2002, Stanway et al. 2003). Opti- 
cal spectra of the brightest quasars at z>6 show a 
clear Gunn & Peterson (1965) effect, demonstrat- 
ing that these objects are located at the end of 
cosmic reionization (Becker et al. 2001, Fan et al. 
2003, White et al. 2003). It is of paramount im- 
portance to study these sources in detail to mea- 
sure the stellar and gaseous constituents, chemi- 
cal abundances and dynamical masses of the host 
galaxies, which in turn constrain galaxy evolu- 
tion models in the very early universe. Studies of 
the host galaxies of bright quasars are extremely 
challenging in the optical and/or near-infrared as 
the central active galactic nucleus (AGN) greatly 
overshines its surrounding. Radio observations of 
the molecular gas phase, however, are a promising 
method to measure such important properties of 
distant quasar host galaxies as the total quantity 
of molecular gas and the dynamical masses. 

The highest redshift quasar currently known is 
SDSS J114816. 64+525150.3 (hereafter: J1148+5251) 
at z=6.42 (Fan et al. 2003). Jl 148+5251 is a 
highly luminous object which is thought to be 
powered by mass accretion onto a supcrmassive 
black hole of mass 1-5 xlO 9 M (which accrets at 
the Eddington limit, Willot et al. 2003). Ther- 
mal emission from warm dust has been detected 
at millimetre wavelengths, implying a far infrared 
(FIR) luminosity of 1.3 x 10 13 L o (Bertoldi et 
al. 2003a), corresponding to about 10% of the 
bolomctric luminosity of the system. J1148+5251 
falls on the radio-FIR correlation found for nearby 
starburst galaxies (Carilli et al. 2004); if the heat- 
ing of the dust is dominated by young stars, 
then the implied star formation rate is an im- 
pressive ~ 3000 Mq yr _1 . High resolution HST 
imaging does not reveal multiple imaging of the 
central AGN (White et al. 2004), rendering strong 
magnification unlikely (cf. the multiple opti- 
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cal imaging in the gravitationally lensed system 
PSS2322+1944, Carilli et al. 2003). 

Using the Very Large Array (VLA) and the 
Plateau de Bure interferometer (PdBI), we have 
detected various rotational lines of carbon monox- 
ide (CO), the most common tracer of molecular 
gas, in Jl 148+5251 (Walter et al. 2003, Bertoldi 
et al. 2003b). This discovery demonstrated that 
abundant, dense and warm molecular gas (the req- 
uisite fuel for star formation) is already present at 
these large lookback times (a redshift of z=6.4 cor- 
responds to an age of the universe of 870 Myr in a 
cosmology with Ho = 71 kms -1 Mpc -1 , £!a=0.73, 
f2 m =0.27 which we adopt in the following). 

Here we present follow-up high-resolution ob- 
servations of the molecular gas distribution of 
J1148+5251 obtained with the NRAO VLA 1 . 
These observations resolve, for the first time, the 
host galaxy of this luminous quasar on kiloparsec 
scales (lkpc = 0.18" at z=6.42) and enable us 
to constrain its molecular gas properties and the 
dynamical mass. 

2. Observations 

Using the VLA in D configuration, Wal- 
ter et al. (2003) detected the CO (3-2) line 
(^ rost = 345.796 GHz) in J1148+5251, which is 
shifted to 46.61GHz (6.43 mm) at z=6.419. The 
emission was unresolved in the D-array obser- 
vations at 1.5" resolution but the line strength 
of 0.6 mJy indicated that follow-up VLA observa- 
tions at higher spatial and spectral resolution were 
possible. In addition, surface brightness argu- 
ments yielded a minimum size of 0.2"(Tb [50K])~ - 5 , 
a resolution which can be reached with the VLA 
B-array. We therefore performed follow-up ob- 
servations of Jl 148+5251 at the VLA in B array 
(60 hours) and C array (24 hours) from January 
- April 2004 (total time on source: ^60 hours). 
Observations were done in fast switching mode us- 
ing the nearby source 11534+49311 for secondary 
amplitude and phase calibration; the phase stabil- 
ity in all runs was excellent (~ 10° phase rms on 
10km baselines). 

Two 25 MHz intermediate frequencies (IFs, cen- 
tered at 46.603250 GHz and 46.622050 GHz, re- 
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spcctivcly) with 7 channels each (channel width: 
3.1250MHz) were observed simultaneously, lead- 
ing to a total coverage of 37.5 MHz (or 240 km s -1 , 
excluding the overlap and edge channel) . This ve- 
locity range encompasses most of the CO line- 
width measured by the VLA and the PdBI 
(FWHM=280kms _1 ), but docs not include the 
line wings and the continuum (in earlier observa- 
tions we derived a 2a upper limit of the continuum 
at 46 GHz of 0.1 mJy, Walter et al. 2003). Given 
the limited signal to noise ratio in individual chan- 
nels, three individual channels were binned to one 
channel, resulting in a total of 4 independent chan- 
nels (width=9.375MHz, or GOkms" 1 ). 

For our analysis we created two datasets. The 
first dataset is based on the new C and B ar- 
ray data with a UV taper of 1 MA, resulting in 
a beamsize of 0.35" x 0.30" and an rms noise 
of Se^ybeam" 1 in a 9.375 MHz channel. The 
rms noise in the map covering the entire bandpass 
(37.5 MHz) is 43^Jy beam -1 . To boost the reso- 
lution in the central region, the B array data only 
were mapped using 'Natural' weighting; this re- 
sults in a beam of 0.17" x 0.13" and an rms noise 
of 45^iJybeam _1 . The original D array data were 
not included in this analysis as they were not ob- 
served in spectral line mode. 



3. Results 

3.1. Global Distribution of Molecular Gas 

In Fig. 1 we present the CO(3-2) emission 
in J1148+5251 over the total measured band- 
width (37.5MHz, 240 kms" 1 ) at 0.35" x 0.30" 
(1.9kpcxl.7kpc) resolution. The emission is 
clearly extended and Gaussian fitting in the map 
plane gives a deconvolved major axis (FWHM) 
of 0.65±0.12" (^3.6 kpc), a marginally resolved 
minor axis of 0.25±0.12" (~1.4 kpc) and a posi- 
tion angle of 15±10° (measured east from north). 
Molecular gas can be seen out to distances of 0.42" 
from the centre (~2.5 kpc, numbers are decon- 
volved for the beam size). The source is possibly 
extended towards the north (see also Fig. 2, first 
channel map) but more sensitive observations are 
needed to confirm this extension. If we assume 
that the main molecular gas concentration forms 
an inclined disk, this gives an inclination of ~65° 



(where 0° corresponds to face-on) which implies 
that any inclination correction to the measured 
rotational velocities is minimal. The fitted peak 
to the central distribution is 0.21 mJybeam -1 
and the integral intensity is 0.55 mJy. This corre- 
sponds to an H 2 mass of M(H 2 ) ~ 1.6 x 10 10 M Q . 2 
A 1 mJy continuum point source 40" north-east 
of the QSO is unresolved, emphasizing that phase 
decoherence is not the reason that the emission of 
J1148+5251 appears extended. 

Given our beamsize of 0.35" x 0.30", the peak 
brightness of 0.21 mJybeam -1 corresponds to a 
surface brightness of 1.1 K, or a beam-smoothed 
brightness temperature of l.lx(l + z)=8.3K at 
z=6.42. 

3.2. Dynamics of the Molecular Gas 

Assuming that the molecular gas is gravita- 
tionally bound and forms an inclined disk with 
a radius of 2.5 kpc (sec previous section) and 
that the gas seen in the PdBI spectrum emerges 
from the same region (full width at zero inten- 
sity: 560 km s -1 or v rot <~280kms , Bcrtoldi 
et al. 2003b), we derive an approximate dynam- 
ical mass for J1148+5251 of - 4.5 x 10 10 M Q (or 
~ 5.5 x 1O 1O M0 if we correct for an inclination 
of ~65°) with an error of order 50%. Within the 
large uncertainties, this number is compatible with 
the derived molecular gas mass and the mass of the 
central black hole, but does not leave much room 
for additional matter (e.g. a massive stellar bulge, 
see discussion below). 

Using the velocity information in our new ob- 
servations, we can now investigate the spatially 
resolved dynamics in the centre of J1148+5251. 
Fig. 2 shows four channel maps (width of one 
channel: 9.375 MHz, or 60kms~ 1 ) which cover 
the entire observed bandpass at 0.35" x 0.30" 
(2.0kpcxl.7kpc) resolution. Due to the intrin- 
sic faintness of the emission, the channel maps 
are noisy, however some interesting features can 
still be identified. In particular there is some in- 
dication that emission is moving from the north 
to the south with increasing channel number (i.e., 
increasing frequency). Observations with better 

2 Here and in the following we assume constant brightness 
in the CO(3— 2) and CO(1-0) transitions and a conversion 
factor to convert CO luminosities to H2 masses of a ~ 0.8 
M (Kkms-'pc 2 )" 1 , see Walter et al. 2003. 
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sensitivity (and covering more of the total CO 
line-width) are required to further constrain the 
internal dynamics of J1148+5251. 

3.3. Distribution of CO at 1 kpc scale 

The distribution of molecular gas at even 
higher resolution reveals that the situation is 
likely more complex than described by the sim- 
ple disk model above. In Fig. 3 we show our 
highest resolution CO (3-2) map derived from 
B-array data only. Here the resolution is only 
0.17" x 0.13", or ~0.95kpcx0.72kpc. 3 This 
map shows that the emission breaks up into 
two regions, a northern source at ll h 48 m 16.640, 
52°51'50.51 (peak flux density: 192+45 £tJy), and 
a southern source at ll h 48 m 16.641, 52°51'50.21 
(flux density: 180±45/iJy), i.e., separated by 
~0.3" (1.7kpc). On a Kelvin scale, both sources 
are at ~4.5K, which corresponds to a lower limit 
for the beam-averaged brightness temperature of 
-35 K at z=6.42. 

To interpret this morphology and its relation to 
the optical quasar, astrometric precision to better 
than 0.05" is needed. However, at the scales under 
consideration it proves to be extremely difficult to 
tie the radio and optical frames together. This 
situation is further complicated by the fact that 
J1148+5251 has a very atypical spectrum (show- 
ing a complete Gunn Peterson effect in the opti- 
cal and only covering parts of the Sloan z band). 
Thus differential chromatic aberration would pre- 
clude tying the SDSS z-band image of Jl 148+5251 
to the radio frame to better than a few tenths of 
an arcsecond even if common sources existed in 
both frames. This uncertainty is indicated by the 
size of the cross in all figures. As the southern CO 
emission is situated more towards the centre of 
the entire molecular gas distribution (cf. Fig. 1) 
we consider it likely, for symmetry reasons, that 
the optical QSO is associated with the southern 
source. In this picture the northern source may 
possibly be the nucleus of an object in the process 
of merging with the quasar host. This situation 
may be reminiscent of the molecular gas distri- 
bution around the z=4.7 QSO BPJ 1202-0725 in 
which two molecular gas peaks have been detected 

3 We note that this is comparable to resolutions typically 
achieved with single-dish telescopes when studying galaxies 
in the nearby universe in the CO(1-0) transition. 



as well (albeit on a factor 10 larger scale, Omont 
ct al. 1996, Carilli et al. 2002). Alternatively, the 
northern source could be a region of high surface 
brightness within J1148+5251. 

4. Summary and Discussion 

We present the first resolved maps of a system 
located at the end of cosmic reionization at z=6.4. 
Our high-resolution VLA data enable us to resolve 
the molecular gas emission in the host galaxy both 
spatially and in velocity space. 

The molecular gas distribution in Jl 148+5251 
is extended out to radii of 2.5 kpc. The central 
region is resolved and shows 2 peaks, separated 
by 1.7 kpc; they account for about half of the to- 
tal emission, with the other half present in the 
more extended molecular gas distribution. Each 
of these peaks harbors a molecular gas mass of 
— 5 x 10 9 M within a radius of 0.5 kpc, respec- 
tively; this mass is similar to the total mass found 
in nearby ULIRGS such as Mrk273 or Arp220 
(Downes & Solomon 1998). The peaks have in- 
trinsic brightnesses of — 35 K (averaged over the 
1 kpc-sized beam) , similar to what is found in 
the centres of nearby active galaxies (Downes & 
Solomon 1998, assuming constant surface bright- 
ness down to the CO (1-0) transition), albeit mea- 
sured over a larger physical area. 

We have assumed that the gas is gravitationally 
bound and is situated in a disk, although the data 
do not rule out the presence of an unbound sys- 
tem, such as an ongoing merger. Based on the ex- 
tent of the molecular gas distribution and the line- 
width measured from the higher CO transitions 
we derive a dynamical mass of — 4.5 x 1O 1O M 
(~ 5.5 x 1O 1O M if we correct for an inclina- 
tion of i — 65°). This dynamical mass estimate 
can account for the detected molecular gas mass 
within this radius but leaves little room for other 
matter. In particular, given a black hole mass 
of mass -1-5 xlO 9 M (Willot et al. 2003), this 
dynamical mass could not accomodate an order 
few xlO 12 M stellar bulge which is predicted by 
the present-day Mbh — Obuigc re l a ti° n (Ferrarese 
& Merritt 2000, Gebhardt 2000), if this relation 
were to hold at these high rcdshifts. Even if 
we assume a scenario in which this bulge was 10 
times the scale length of the molecular gas emis- 
sion, we would still expect a bulge contribution of 
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fewxlO 11 M within the central 2.5 kpc (assum- 
ing a central density-profile of p ~ r~ 2 , e.g., Jaffc 
1983, Tremaine 1994) which can not be reconciled 
with our results. Our finding therefore suggests 
that black holes may assemble before the stellar 
bulges. This would be in conflict with the popu- 
lar picture of co-eval evolution of the central black 
hole and the stellar bulge. The derived dynami- 
cal mass and the hypothetical stellar bulge mass 
(based on Mbh — "bulge) can be only brought in 
agreement if (i) our derived inclination is signifi- 
cantly off (i.e., i < 10°), (ii) the estimated black 
hole mass is over-estimated by orders of magni- 
tudes (which would imply highly super-Eddington 
accretion) , or (iii) if the assumption that the gas is 
gravitationally bound is seriously wrong. We note, 
however, that depending on the space density of 
similar sources at these high redshifts, the smaller 
dynamical masses suggested in this study may be 
in better agreement with the masses predicted by 
CDM simulations in the very early universe as a 
10 12 Mq bulge would imply a rather massive dark 
matter halo of > 1O 14 M . 

This study shows the potential for future stud- 
ies of the molecular gas content and dynami- 
cal masses in the highest redshift galaxies us- 
ing ALMA, where resolutions of < 0.1" will be 
achieved routinely. 

We would like to thank the anonymous referee 
for very helpful comments. 
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Fig. 1.— CO(3-2) map of J1148+5251 of the 
combined B- and C-array dataset (covering the 
total bandwidth, 37.5MHz or 240kms~ 1 ). Con- 
tours are shown at -2, -1.4, 1.4, 2, 2.8, 4xcr 
(lo-^^Jybeam- 1 ). The beamsize (0.35 x0.30") 
is shown in the bottom left corner; the cross indi- 
cates the SDSS position (and positional accuracy) 
of J1148+5251. 
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Fig. 2. — Channel maps of the CO(3-2) emis- 
sion based on the combined B and C array data 
sets. The same region is shown as in Fig. 1 (Dec- 
lination is in arcsecond relative to 52° 51' and 
the Right Ascension, labeled on the third chan- 
nel, is relative to llh 48m). One channel width 
is 9.375 MHz, or 60kms _1 (frequencies increase 
with channel number and are shown at 46.597000, 
46.606375, 46.615750 and 46.625125 GHz). Con- 
tours are shown at -2.8,s -2, -1.4, 1.4, 2, 2.8, 4xcr 
( 1 er=86^ Jy beam" !). The beamsize (0.35 x 0.30") 
is shown in the bottom left corner; the cross indi- 
cates the SDSS position (and positional accuracy) 
of J1148+5251. 
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Fig. 3. — CO(3-2) map at ~lkpc resolution of 
the central region displayed in Figs. 1 and 2. Con- 
tours arc shown at -2.8, -2, -1.4, 1.4, 2, 2.8, 4xcr 
((7^45/iJybeam- 1 ). The beamsize (0.17 x 0.13") 
is shown in the bottom left corner. The positional 
uncertainty of the SDSS position is of order ±0.1" 
or about the size of the cross. 
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